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^ ■ We prove the generalized induction equation and the generahzed local induction equa- 

pL| , tion (GLIE), which replaces the commonly used local induction approximation (LIA) 

to simulate the dynamics of vortex lines and thus superfluid turbulence. We show 
that the LIA is, without in fact any approximation at all, a general feature of the 
Ph! velocity field induced by any length of a curved vortex filament. Specifically, the LIA 

^— » ' states that the velocity field induced by a curved vortex filament is asymmetric in the 

binormal direction. Up to a potential term, the induced incompressible field is given 
by the Biot-Savart integral, where we recall that there is a direct analogy between 
^ ■ hydrodjTiamics and magnetostatics. Series approximations to the Biot-Savart inte- 
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grand indicate a logarithmic divergence of the local field in the binormal direction. 
While this is qualitatively correct, LIA lacks metrics quantifying its small parame- 



(N 

O ' ters. Regardless, LIA is used in vortex filament methods simulating the self-induced 

motion of quantized vortices. With numerics in mind, we represent the binormal 
field in terms of incomplete elliptic integrals, which is valid for M^. From this and 

X ■ 

^ . known expansions we derive the GLIE, asymptotic for local field points. Like the 

LIA, generalized induction shows a persistent binormal deviation in the local-field 
but unlike the LIA, the GLIE provides bounds on the truncated remainder. As an 
application, we adapt formulae from vortex filament methods to the GLIE for future 
use in these methods. Other examples we consider include vortex rings, relevant for 
both superfluid ^He and Bose-Einstein condensates. 
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I. INTRODUCTION 

The term superfluid denotes a phase of matter whose dynamical flows can be described, 
at flnite non-zero temperature, by a two-component macroscopic fleld with well-deflned 
properties." One component is a purely classical fleld, while what remains is called the 
superfluid component. The superfluid is ideal in the sense that it is inviscid and has inflnite 
heat capacity provided by its lack of classical entropy. Rotation enters the superfluid compo- 
nent in quantized vortex fllaments^i^ that, for example in ^He, transmit thermal information 
acoustically and are detected by this second-sound.- Superfluid dynamics can be generated 
by the introduction of a small heat flux.- Conservation of mass requires that the classical 
movement away from the heat flux be offset by a counterflow of the superfluid component. 
If this heat flux is not small, then the quantized vortices tangle, indicating the onset of 
superfluid turbulence.-""— For large heat fluxes, the superfluid transitions into a purely clas- 
sical phase. In classical turbulence, vorticity can concentrate into complicated geometries. 
For this reason, large-scale simulation of classical vortex dominated flows is computationally 
costly. 

Vortex line structures are most appropriate to superfluid models of ^He where the quan- 
tized filaments have radii of a few angstroms.-^^^ These quantized vortices provide a coherent 
structure for aggressive analytical and numerical study unavailable to classical fluids. Vor- 
tex line structures can also be used to model atomic Bose-Einstein condensates, where there 
has recently been a revival of interest in superfluid turbulence due in part to a series of 
remarkable experiments in the Bagnato group.— Our study begins with simpliflcations to 
the Biot-Savart representation of the fleld induced by a vortex line 



The associated reduction of dimension aids analytic calculation and reduces numerical cost. 
When such fllaments are considered initial-data to the Navier-Stokes problem, then global 
well-posedness results.— The use of these data to approximate self-induced vortex motion 
is the backbone of the vortex filament methodM^ The question of numerical convergence 
and accuracy of such vortex fllament techniques has been addressed affirmatively in the 
literature.— Vortex filament methods reduce the cost of large-scale simulations by restricting 
analysis to the local field. Due to the complexity of classical vortical fiows, interest in these 
methods waned during the 1980s.— '^^i'^"— However, filament methods are highly appropriate 
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for the constrained vortex structure associated with a superfluid. Although ([T]) provides a 
straightforward starting point for numerical computations, vortex filament methods avoid 
numerical integration altogether by replacing ([T]) with a local induction approximation. 

The local induction approximation (LIA) is a result from classical fluid dynamics, which 
states that a space-curve vortex defect of an incompressible fluid field with nontrivial cur- 
vature generates a binormal asymmetry in the local velocity field. That is, the field local to 
a length of curved vortex filament induces a flow, which generates fllament dynamics. This 
result, known by Tullio Levi-Civita and his student Luigi Sante Da Rios in the early 1900's,— 
was rediscovered by various post World- War II groups.—""— Together, Ricca^^ and Hama^ 
provide an excellent chronology of LIA, a topic now common in vortex dynamics texts.— i^*^ 
Exploration of LIA occurs in various settings including differential geometry,—"— differential 
equations2^ and limits of matched asymptotic expansions of vortex tubes.— i^^ii^ 

Our derivation avoids the complications of matched asymptotic expansions by treating 
vorticity concentrated to an arc. Under this geometry Eq. ([1]) can be reduced to a canonical 
elliptic representation without the use of power-series approximation to the Biot-Savart 
integrand.—!^ While Taylor approximation quickly reveals binormal flow as a dominant 
feature of the induced fleld, it lacks error bounds and is often restricted to a two-dimensional 
subspace of M^. This paper resolves both issues by recasting Eq. ([1]), for a vortex-arc, 
into an elliptic form valid for all field points. Using this form, one can then use a known 
asymptotic formula to represent the local field. We offer that this should be adopted, instead 
of LIA, for use in vortex filament methods and Schwarz's description of the Magnus force.— 
Specifically, we will prove the following results: 

Theorem 1 Generalized Induction Equation 

Let u: = V X V be localized to an arbitrary arc with parameterization ^ = {Rsm{6),R — 
i?cos(^),0), where R G and 9 G Dl = {—L,L] for some L G (— 7r,7r]. Then there exists 
bounded functions Vi, a\, a2, L± andk, of s = \x\/R = k,\x\, such that the induced velocity 
field is given by 



where F is an incomplete elliptic integral of the first kind. Moreover, there exist constants 



v{x) = Vi{e) aiie) [F{L+, k) - F{L_, k)] + a2{e) 



dF{L+,k) dF{L_,k) 



de de 
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Pi, ^2, Ps, Pa such that Vi can be written as 

V^{e) = ePit - eP2h+ {eP^ + eP^ + P^) b (3) 
where t, h, b, are the tangent, normal and binormal vectors of the local coordinate system. 
Theorem 2 Generalized Local Induction Equation 

Under the same hypotheses of theoremUl and for e -C the induced velocity field is dominated 
by the binormal flow, 

vs{x) = Anx, [a,{e) [F{L^, k) - F{L_, k)] + a,{e) ^Ei^llR _ ^^^^f^ ) b. (4) 

where X2 is a dimensionless angular component of the spherical decomposition of x and k 
is the curvature of the vortex arc. The limits e — and L — )■ imply that A; — j- 1 and 
A — 7- and in this case the incomplete elliptic integral of the first kind admits the asymptotic 
relation F ~ Fi where 



Using this, along with standard differentiation formula for incomplete elliptic integrals of the 
first kind, provides a first order asymptotic form for the local field given by 



9x2^1 (A, k) X2E{L, k) ^ sin(2L) \ ^Jl + k^ sm\L) + ^/T^W^^^ 



2 {l-k^)k 2{l-k^)^/l-k^sm\L) 

(6) 

where E is an incomplete elliptic integral of the second kind. 

In words, the first theorem expresses the velocity field generated by vortex arc in terms 
of incomplete elliptic integrals of the first kind. Moreover, this field can be decomposed 
into three fields controlling the tangential, circulatory and binormal fiows. Of these fields, 
the binormal contribution is 0(1) while the remaining fields are 0{e). The second theorem 
considers the remaining field in the limits of e — ?■ and L — j- 0. In this limit the incomplete 
elliptic integral of the first kind admits an asymptotic form and consequently provides a 
representation for the velocity field local to the vortex arc. This asymptotic form is com- 
parable to LIA in that the Biot-Savart integral has been 'resolved' and the binormal fiow is 
represented by elementary functions. This form is only valid for filaments of infinitesimal 



arclength and consequently idealized. However, using the same asymptotic framework, one 
can construct expansions valid for arcs of finite length. In fact, the remainder terms of such 
expansions are known and thus the associated approximation error can be controlled. The 
necessary asymptotic results are quoted in the appendix. 

The rest of this document will be organized as follows. In Section II we define the 
geometry and derive the Biot-Savart representation of the induced velocity field. In Section 
III we convert this representation into an elliptic form and prove the generalized induction 
result ([2])-(|3]). In Section IV we reduce this elliptic form into a sum of incomplete elliptic 
integrals of the first kind. Lastly, using known asymptotic results, we derive an expression 
for the local velocity field and prove the generalized local induction equation (GLIE) result 
We conclude with some discussion on adapting this result to vortex filament methods 
and prospective avenues of future work. 

II. BIOT-SAVART AND QUANTIZED VORTEX RINGS 

It is well known that a vortex-defect with trivial curvature embedded into an incompress- 
ible fluid does not induce autonomous dynamics. This is due to an angular symmetry in 
the induced velocity field. This symmetry is no longer available for curved vortex elements. 
Using a vortex ring, it is possible to introduce nontrivial curvature and avoid approximations 
to the Biot-Savart integral. To be precise, we treat a vortex structure a; : — t- M'^ such 
that 



for K ^ = R E Dl = {—L,L] and L G [0,7r]. Thus, at the point x = {xi,X2,X3) we get 
an element level description of the velocity field. 




(7) 



where ^ : D — )■ M^, C M, is parameterized by the ring 




(8) 
(9) 



d^= Rcos{9)i + Rsm{9)} d9 




(10) 



D, [1x12 + 1^12-2(51^1 + 526+536)] 
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where we have used the Levi-Civita symbol, eijk, and employed silent-summation over re- 
peated indices. Noting that = 2R^ — 2R'^ cos{6) provides the formulae 

, , 9 f sin(6') , , 

v^ = -\^\k\s -^de, (11) 

V2 = \^\^'x,l^ "-^de, (12) 

n, , f xi sin(6') — X2 cos(^) ,^ f cos(6') — 1 ,^ 

= " X, ^^37^ + " I, -^D^"" (l^*) 

where the denominator is given by 

D = [ci + C2 cos(0) + C3 sin(e)] (14) 

R^ci = |x|2 + 2R^ - 2x2/?, (15) 

i?2c2 = 2|x|x2i?- 2i?2, (16) 

R^cs = -2\x\xiR. (17) 

For future limiting work, we have chosen a radial-representation for x = |x|(xi, X2, X3) where 
is the i^^ dimensionless angular component of x. 



and whose coefficients are 



X 



III. CONVERSION TO ELLIPTIC FORM 

The previous integral representations for the velocity field can be cast into elliptic form. 
To do this, we first reduce each integral into an elliptic integral by taking derivatives with 
respect to internal parameters. Doing so gives 

^ d f d9 

vi = 2KXse— / — =, (18) 
dc3 Jdl vD 

d r de 



V2 = -2KX3e— / -j=, (19) 
dc2 Jd, VD 



V3 



d „ d d ^ d 

2Kex2- 2nexi- h 2n- 2n- — 

dc2 dcs dci dc2 



(20) 



where the parameter e = |x|/i? is the ratio of radial distance to the radius of curvature. 
Application of the chain-rule gives the induced velocity field as 

/X ^ ^ , s. d f d6 , . 
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where the vector Vi is given by 

(22) 



Vi(£:) = 2eKX3—-^i — 2eKX2—r-] 

de de 



dc2 dcs dci dc2 

2Kex2—, 2Kexi— — h 2k— 2k—— 

de de de de 



implying that the velocity field is determined by the derivative of an incomplete elliptic 
integral. Moreover, this proves Eq. (jS]) from theorem [1] where 

Pi = 2«:X3^, (23) 

rlr 

[32 = 2«:a;2^, (24) 

dc'^ , . 

= -2«:xi^, (25) 

For e ^ 1 we find that the velocity field is dominated by 

/ N „ d f dO s- 
..(x)^-8...-y^^^k, (27) 

Figure [1] shows the vortex configuration as well as the associated field and vortex coordinate 
geometry. Using the depicted spherical decomposition of x we find that the previous dimen- 
sionless parameter is given by X2 = sin(7i) sin(72). Moreover, we observe that the standard 
basis vector k corresponds to the binormal vector b. These two facts show that the velocity 
field, asymptotically close to the vortex arc, is asymmetric in the binormal direction and 
that this affect is extremized for field-points on the normal-axis, which agrees with standard 
results of induced binormal flow. 



IV. REDUCTION OF ELLIPTIC FORM TO CANONICAL ELLIPTIC 
INTEGRALS 

Before we construct the asymptotic representation of the velocity fleld, the previous 
integrals are converted into canonical forms. The induced velocity fleld is controlled by an 
integral of the form 

de 



(28) 

+ C2 cos(6') + C3 sin(6') 

which can be converted to a sum of incomplete integrals of the flrst kind. First, we introduce 
a new angle deflned by tan(0) = C3/C2 and hypotenuse r"^ = cl + c| to get 

de de 



Dl ^/ci + C2 cos{e) + C3 sin(^) J-L ^/ci + r cos 
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(a) Vortex Configuration (b) Tangent-Normal Plane 

FIG. 1: Global and Local Coordinate Geometry 
In subfigure (a) the vortex arc is depicted in M'^ where the circle parameterization, C, is 
composed of the solid line representing the vortex filament and a dashed line representing 
a continuation of the parameterization. These two regions are separated by the cut-off 
parameter L. This subfigure also shows the spherical decomposition of the field point x 
where 71 is the azimuthal angle and 72 is the polar angle associated with the spherical 
decomposition of x. Lastly, this subfigure shows the configuration of the Serret-Frenet 
local basis vectors t, fi, b, which, for ease of use, are oriented to correspond to the standard 
global basis vectors for M^. In subfigure (b) the projection of subfigure (a) onto the x-y 
plane is given and shows the polar decomposition of the filament point 
$, = {Rsm{9),R- Rcos{9),0). 



Now, introducing a change of variable 2tlj = (j) — 6 and the notation L± = ((/) ± L)/2 we 
apply trigonometric formulae to get 

r = ^ r (30) 

where /c^ = 2r/(ci + r). Lastly, 

de _2[F{L+,k) - F{L_,k)] 

Dl v^ci + C2 cos(^) + C3 sin(^) y/cTTr 



(31) 
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where F is the standard incomplete elhptic integral of the first kind, 

F(j,, k) = r ■'^ = r (32) 



such that A = sin((y9). 



V. ASYMPTOTICS FOR THE INCOMPLETE ELLIPTIC INTEGRAL OF 
THE FIRST KIND 

Having reduced the Biot-Savart representation of the velocity field to a canonical form, 
we can now make use of the known asymptotic formula of Karp and Sitnik,— which permits 
the study of f l32|) for all (A, k) G [0, 1] x [0, 1]. Specifically, they derive a series representation 
and remainder term for F, which is asymptotic for — > 1. Their complete theorem is quoted 
in the appendix, but we only require the first-order approximation 



which is asymptotic to F for A —t- and k 1. This asymptotic formula is not suited to 
differentiation.—"— Thus, we must first apply the differentiation formula, for the incomplete 
elliptic integral of the first kind, prior to its asymptotic evaluation. Doing so gives lengthy 
formulae and for these we introduce the constants 

A = ^I^IH^, (34) 
r 

A, = j-^^^{e-X2 + A), (35) 

A2 = ^i=, (36) 

A3 = -2^^-^^, (37) 



A. - + ( ^^^^f'-\f^ ) A. (38) 

(ci+r)3/2 VrV (ci + r)2 J ^ ^ 

Using these constants, find 

^ [ 4^ = 2Ai[F(L+, k) - F{L., k)] + A2 [fi(A3, A4, As, k, L+) - 1](A3, A4, As, k, L_)] 
de Jdl VD 

(39) 
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where 

A3 A^E{L,k) AiF{L,k) Ai2ksm{2L) 



n{As,A4,A5,k,L) 



Vl-P sin^(L) 1(1 -k^)k 2k 4{l-k^)^l-k^ sm\L) 

(40) 

is given by differentiation formula for incomplete elliptic integrals of the first kind. Together 
with Eq. (pTl) . proves Eq. (j2]) of our first theorem where ai = 2Ai and 02 = ^2- From Eq. 
(l27|l we find that the local velocity field is given by 

v,{^) = -8KX2 (2Ai[F(L+, k) - F{L_, k)] + A2 [n{A3, A4, As, k, L+) - n{As, A4, A5, k, L_)]) . 

(41) 

At this point, the asymptotic formula (lA.ip can now be applied to F{X, k) where A = 

sin(L-i-). To compare these results to standard LIA we take e — and L — 0. In this case 

Ci = — C2 = r ~ 2, C3 ~ and the constants take the asymptotic forms 
X2C2-X1C3 

A = ~ -X2, (42) 

r 

^1 = ( 7^^3/2 - X2 + A) ~ 0:2, (43) 

(ci + r)'^/^ 

A2 = ~ 1, (44) 

. ^ X2C3 + XiC2 . . 

A3 = -2 -3 xi, (45) 

= (c, + r)3/2 + V - ^ ^T^Tp J ^ ~ - T- 

Together this gives the first-order asymptotic representation for the velocity field, 

9x2Fi(A, k) X2E{L, k) ^ si^(2^) [V^ + sin^(^) + v^l-Zc^ sin^(L) 



VelXl ~ —6KX2 



2 (1-A;2)A; 2(1 - P) v^l - A;^sin4(i:) 

(47) 

for the limits e — 0, A; — 1 and L — )■ 0, A — ?■ 0. This proves Eq. ([6]) of our GLIE. It should 
be noted that the above formula is nonzero even for the extreme case ofL = l — A;— J-O. 
The physical meaning of this statement is that the local field induced by an infinitesimal 
segment of a vortex line is nonzero and asymmetric in the binormal direction. 

VI. DISCUSSION AND CONCLUSIONS 

We have derived an asymptotic representation for the local velocity field induced by 
a curved vortex filament. This derivation generalizes the previously known statements of 
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induced binormal flow, which play an important role in two-component superfluid simulation. 
In such simulations one must calculate the superfluid and normal fluid flows as well as 
their mutual friction interaction. This mutual friction embodies the scattering of rotons 
and phonons off of the vortex structures.-"—'^ It is possible to calculate this interaction 
in a manner self-consistent with Navier-Stokes and fully coupled to both components.— 
The normal fluid is approximated through Navier-Stokes simulation techniques while the 
kinematics of the superfluid make use of LIA. Though our focus is LIA dynamics we offer 
the following references to the computational ffuid dynamics literature, which has been 
used for coupled two-component superfluid simulations.—"— While there has been progress 
in these techniques,—!^ the recent growth of the highly adaptable discontinuous Galerkin 
methods^ and their application to nonlinear fluid flow and acoustic problems^ is especially 
provocative. 

Mathematically, the kinematics of a vortex filament, ^, are described by^ 



^ = V5 + V, + /3^'x (Vjv-V5-V,)-/3'|'x [I'x (V^-V5-V,)]. (48) 



In a filament method, it is typical to prescribe the normal-fluid flow Vat and neglect the 
mutual friction terms involving /3 and /3'.— This leaves only a potential flow V5 and induced 
flow V/. Of the remaining quantities, the computationally costly induced flow is managed 
through the LIA, 



where L± = (0±L)/2 is related to the cutoff length L and angle </>. Vortex fllament methods 
avoid integration by application of this approximation to nodal points of the Lagrangian 
computational mesh attached to the fllament centerline. Alternatively, we could simply 
replace LIA with GLIE dl])-© and write V/ ^ Vg and prescribe a fleld point x and arclength 
s = 2RL. However, if the higher-order circulatory and binormal terms are desired, one could 
use ([S])-© and employ efficient numerical routines for the incomplete elliptic integrals.—"— 
Either of these changes will then be applied to node points of a computational mesh modeling 
the filament structure. 

The use of piecewise linear interpolants, while prevalent in numerics, cause spurious effects 
when applied to a vortex centerline. The interpolants themselves have zero local curvature 
and their connections form cusps with undefined local curvature. Typically, local induction is 




(49) 
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applied to higher-order interpolations. While one can use the generalized induction equation 
or GLIE on this mesh, the natural vortex-arc construct has been adapted to efficient meshing 
techniques.— Consequently, computational cusps are avoided and local curvature is always 
well-defined when GLIE is applied to such vortex-arc meshes. Lastly, what remains is re- 
meshing to allow for the experimentally witnessed vortex nucleation.— 

Meshing is the most difficult aspect of vortex filament implementations. Not only must 
the mesh adapt to the vortex dynamics, it must be made to reconnect filament elements that 
are not predicted by the Eurelian theory.— The most elementary reconnection algorithms 
appeal to nonlinear Schodinger theory and force reconnection of filaments passing within 
a few core widths of each other.— The current theory of the reconnection process is not 
satisfactory and efforts to avoid ad hoc simulated reconnection continue.—"— 

Superfluid turbulence dominated by quantized vortex flows is an active area of analytic, 
numerical and experimental research.— i^"— Though local induction techniques will play a 
part in continued numerical investigations, understanding geometric and topological quan- 
tiflcation of a tangled state is as important and still a work in progress.—"— Lastly, the 
vortex line approximation, while useful and appropriate, must eventually be discarded in 
favor of nontrivial core-structure. It is likely that the methods developed within this paper 
can be adapted to current arguments used to study flelds induced by vortex tubes.— 

That being said, this work makes it clear that binormal flow proportional to curvature is 
a general feature of vortex fllament dynamics. This means that the well-celebrated trans- 
formation of Hasimoto^, which connects the fllament 's curvature and torsion variables to a 
wavefunction controlled by nonlinear Schrodinger evolution, is fundamental to vortex flla- 
ment dynamics. Consequently, even geometrically complicated fllament dynamics are rooted 
in integrable systems theory. This connection underpins efforts to predict allowed fllament 
geometries from the associated integrable systems.—"— 
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Appendix: Asymptotic Representation for Incomplete Elliptic Integrals of the 
First Kind 



The following theorem is one of the two major results proven in Karp and Sitnik.— The 
second result gives a simpler expression but is not valid on the leftmost edge of the unit 
square and therefore not used in our calculations. 

Theorem 3 For all (A, k) G [0, 1] x [0, 1] and an integer N > 1, the previous elliptic integral 
admits the representation 



(A.l) 



where is given by the recurrence formulae 

_ anSn+2{x), +bnSn+l{x) + CnSnjx) + kn 

" 4(n + 3)2 ' ^^-^^ 

a„(x) = 8n^ + 36n + 42-x(2n + 5)^ (A.3) 

fe„(x) = 2x(4n2 + Un + 13) - (2n + 3)^, (A.4) 

Cn{x) = -4x(ra + 1)^ (A.5) 
h (r) - ^(2^ + 5)(2r2 + 4)^ + {n + 3){8n' + 24n + 17) 

.o(x) = -21nfi±^V (A.7) 



.(.)^g-l)ln(l±^)-lv^.l.|, (A.S) 

, , / 9 2 X 3\ fl + VTT^\ f9 7\ r- 7 1 21 2 

= ^-32^ + 4 - 4 J [ 2 J + ^32" - 16 J +16 + 8^-64^ 

(A.9) 

and the remainder term is negative and satisfies, 

^pr^yyip k) < -Rn[\ k) < 2[(Ar+ 1)!]2 ^^■^^> 

where the positive function 



In 



v/l+(l-A2)/[aA2(l-fc2)]+l 



fN{\ k) = — <^ , ,^ „ == - (1 - k') In 



1 - a(l - P) aA^l + (1 - A2)/[aA2(l - A;^)] ' V^"^ 



|a=(Ar+l/2)2/(Af+l)2 

(All) 
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is bounded on every subset of E of the unit square, where 



sup 



1 - k 
1 - A 



< oo 



(A.12) 



and is monotonically decreasing in N . 



From this theorem we denote its first-order approximation as 



Fi(A, k) = In 




) 



+ 



1-P 
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In 



( 



1 + A 
1 - A 




and note that this expression is asymptotic in the A variable. 
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